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Radio emission is insensitive to dust obscuration, and the breadth of the radio luminosity function ensures 
that sources are detected over a wide range of redshifts at all radio flux densities. As a result, radio continuum 
observations are an efficient and unbiased probe of both nuclear (AGN) and star-forming activity over all cosmic 
epochs. 

The SKA's "ultra-deep radio continuum surveys" will provide the answers to at least three key astrophysical 
questions: (i) the star-formation history of the Universe, (ii) the evolution of the low-power end of the radio 
galaxy luminosity function and (iii) the relationship between the radio-loud AGN, star-formation and radio-quiet 
AGN phenomena. 

In this paper we discuss the AGN science that will be enabled by the deep radio continuum studies using 
the SKA. It is important to recognise that the sub-/ijy SKA sky will be dominated by populations other than 
'radio-loud' AGN. In this way the SKA will not only provide an unbiased tracer of the star-formation history 
of the Universe but also be able to study the populations of sources we currently describe as 'radio-quiet'. To 
illustrate this point we present simulations of the extragalactic radio sky based from models of the evolution of 
the radio luminosity function. From these simulations we predict typical source distributions and estimate the 
natural confusion limit to ultra-faint flux density limits relevant to the science and design goals of the SKA. 



1. INTRODUCTION 

As a consequence of the broad radio luminos- 
ity function (RLF) which spans many decades 
in radio power, radio continuum observations at 
all flux density limits detect extragalactic sources 
across a wide redshift range. Current large-area 
radio surveys, e.g. the NRAO VLA Sky Survey 
(NVSS, El), Faint Images of the Radio Sky at 
twenty cm (FIRST, |57j) and the Sydney Uni- 
versity Molonglo Sky Survey (SUMSS, 0), cover 
thousands of square degrees of sky to mJy-flux 
density limits. These surveys are the deepest 
practicable over large areas of sky given the large 
amounts of observing time involved. Surveys to 
increasingly fainter flux density limits cover far 
smaller areas, e.g. the Australia Telescope ESO 
Slice Project (ATESP, |2T]) covers 26 square de- 
grees at 1.4 GHz to 5*1.4 GHz = 0.4 mJy (5cr 
rms) whereas the VLA observations of the Hubble 
Deep Field (north) area covers just 0.35 square 
degrees [23 to 61.4 GHz = 0.04 mJy (5er rms). 



The SKA, with its significantly larger collecting 
areas, will be capable of large-area sub-/iJy sur- 
veys as well as deep nano-Jy small-area surveys. 
The challenge is to determine what the sky might 
look like at these faint limits. Firstly it is obvious 
that the SKA will be able to compile a complete 
census of AGN throughout the Universe, map- 
ping their evolution and reveal the relationship 
between the various populations. In turn this will 
allow us to unravel the underlying physics of the 
AGN lifecycle (Jarvis & Rawlings, Falcke et al, 
this volume). The key AGN questions are 

- When do the first AGN form? What is their 
role in the epoch of reionization? 

- What drives the evolution of AGN? Is it in- 
fall of material via mergers, internally-changing 
accretion rates? 

- What is the typical lifecycle of an AGN? Is it 
the same galaxies re-activating? 

- How many populations of AGN are there? 

- Is there a continuum of AGN activity, from 
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the super-Eddington accretors to normal galax- 
ies (and perhaps to even lower-power/low mass 
objects)? Are there objects with BH's of a few 
hundred solar masses? 

- What is the interplay (if any) between AGN 
activity and star formation? 

In this contribution we concentrate on how 
SKA deep continuum surveys will revolutionise 
our study of the AGN phenomenon in light of 
our current knowledge. 

2. RADIO-LOUD AGN, CONTINUUM 
SURVEYS AND THE RLF 

Bright continuum radio surveys reveal sources 
which are powered by non-thermal (synchrotron) 
radio emission from AGN embedded within indi- 
vidual galaxies span a huge range in luminosities. 
AGN activity is detected across a huge radio lu- 
minosity range from between 10 22 W Hz -1 to 10 28 
W Hz" 1 at 1.4 GHz. 

At low radio powers we find classes of objects 
more usually termed radio-quiet AGN: These in- 
clude Seyfert galaxies and QSOs. Because of the 
huge range of observed luminosities, it can be ar- 
gued that a continuum of nuclear radio activity 
exists - such that all galaxies may harbour a cen- 
tral MBH and it is only periodical triggering of 
activity which gives rise to radio jets and lobes 

Radio sources are classified by types, e.g., nor- 
mal galaxies, Seyferts, FRI or FRII radio galaxies 
or quasars, based on their luminosity, radio mor- 
phology, radio spectral index, variability, optical 
counterpart, and IR, optical and X-ray spectral 
features. 

The existence of a continuum of radio activ- 
ity is clear, with the various source populations 
overlapping at all but the highest radio lumi- 
nosities. This is clearly seen in recent, accu- 
rate determinations of the local radio luminosity 
function (LRLF, 7 24 ) for moderate power radio 
sources (10 22 to 10 24 W Hz^ 1 derived from com- 
prehensive galaxy redshift surveys (UGC, 2dF- 
GRS, SDSS) matched to large-area radio sur- 
veys (NVSS, FIRST, SUMSS). At lower radio 
luminosities, the LRLF can only be compiled 
by observing individual objects for long inte- 
gration times. At the high luminosity end the 



space density of objects is small at all epochs, so 
that all that can be determined is the sparsely- 
sampled (evolving) RLF from large-area samples 
(e.g. |2H1 )■ Moreover the high luminosity end of 
the LRLF can only be derived from evolutionary 
fits to the source count, as it cannot be directly 
determined due to the very low local space den- 
sity of these objects. 

In terms of source types, powerful extragalactic 
radio sources are a mixture of radio- loud galaxies, 
quasars and BL Lac objects - all members of the 
Active Galactic Nuclei (AGN) class of objects. 
At the zeroth level, radio galaxies can be classi- 
fied as either Fanaroff-Riley class I or II (FRI or 
FRII): These are distinguished by the distribu- 
tion of bright knots of radio emission and colli- 
mation features ■ In general, FRIs are of lower 
radio power than FRIIs, although there is consid- 
erable overlap in the radio luminosity function of 
the two classes. The current paradigm for radio- 
loud AGN attributes FRI and FRII radio galax- 
ies as the 'parent' populations of the great ma- 
jority of the observed powerful AGN (e.g. |25j . 
such that quasars, blazars and BL Lac objects 
are manifestations of FRI and FRII radio galax- 
ies aligned close to our line of sight, with their 
core radio emission Doppler-beamed (\£jg\ and ref- 
erences therein). 

Typically, FRI and FRII radio sources have lin- 
ear dimensions of the order of between one and 
a few hundred kilo-parsecs. Their radio struc- 
ture shows a pair of oppositely directed radio 
lobes which are sometimes joined to the radio core 
via a pair of collimated, knotty jets of bright ra- 
dio emission. The radio emission from the lobes 
of FRI and FRII radio galaxies is usually op- 
tically thin with a radio spectral index around 
—0.7, although these sources may also have a 
self-absorbed flat-spectrum core. Compact radio 
quasars, and radio galaxy cores, typically have 
angular dimensions much less than an arc-second, 
and due to self-absorption have flat radio spectra. 
The compact radio sources (quasars and blazars) 
are often variable across a wide waveband (radio- 
to-UV) and their host galaxies may show broad 
optical/UV emission line spectra. 

The study of the evolution of the most powerful 
radio-loud AGN is based on the RLF. It is well es- 
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tablished that powerful evolution is required, for 
both these objects, and radio-quiet QSOs: The 
form of this evolution - whether luminosity evo- 
lution 0] or luminosity-density evolution [S] is 
not yet determined, but pure luminosity evolution 
infers source lifetimes which are unphysical |17| . 
Observations at SKA sensitivities are required to 
detect low power radio sources at high redshift; 
this will determine whether all radio-loud AGN 
co-evolved or if we are seeing multiple popula- 
tions undergoing bursts of AGN activity. 

Classification of individual radio sources is of- 
ten hampered by the difficulty of optical/IR/UV 
follow-up: if the optical host is identifiable it can 
still be too faint for 8-m class spectroscopy. For- 
tunately cross-waveband detection may not be a 
requirement for the SKA as the distinction be- 
tween AGN and/or star- formation activity is best 
made via VLBI-scale imaging of /iJy sources. 

2.1. Radio sources and the SKA 

Significant progress has been made in under- 
standing the general radio AGN phenomenon 
from both 

(i) wide and shallow sky surveys, cataloguing up 

to 100 sources per square degree to about 
1 mJy (e.g. FIRST), and 

(ii) small and deep pencil beam surveys (<<1 
degree areas), to very faint flux densities 
e.g. the Hubble Deep Field North and other 
fields ^01 and the Hubble Deep Field South 

The primary gain of the SKA will be its 
ability to reach deep sensitivity limits and 
probe the enormous span of AGN radio powers 
across a wide redshift range. This will remove 
the luminosity-rcdshift (P-z) degeneracy which 
plagues flux density limited samples of today: at 
significant redshifts (z>l) only the highest-power 
sources are detected £Q. 

Multi- wavelength (optical/UV) imaging and 
spectroscopy of sources in the wide and shal- 
low surveys have revealed that the vast major- 
ity of radio sources stronger than a few mJy at 
1.4 GHz are powerful radio galaxies or quasars. 
Moreover, due to the steepness of the high ra- 



dio power RLF, coupled with strong cosmic evo- 
lution between z=l and z=5, the entire popula- 
tion of high power radio AGN is sampled above a 
few mJy: below this flux density limit the nature 
of the radio source population changes dramat- 
ically. Small, deep pencil beam surveys reveal 
that below about 1 mJy at 1.4 GHz, the differ- 
ential radio source count flattens: both the star 
forming galaxy population and galaxies with low 
power AGN contribute. The contribution from, 
and physical properties of the star forming galax- 
ies is discussed in van der Hulst et al. (this vol- 
ume). 

3. THE SKA EXTRAGALACTIC CON- 
TINUUM SKY 

The SKA will offer significant improvement in 
sensitivity so that it will be possible to detect 
radio emission as weak as 10 nJy within ~ 1000 
hours integration time. In this way, the SKA will 
probe a huge part of the very broad RLF at all 
epochs and the misnomers of 'radio-loud' - 'radio- 
quiet' may be buried forever. 

High resolution, multi-wavelength observations 
will be critical to distinguish those extragalac- 
tic radio sources with AGN emission processes. 
Whilst nothing is known about the radio source 
count below about 30/iJy, some estimates can be 
made by extrapolating the RLF to low luminosi- 
ties (10 18 W Hz -1 ) and adopting evolution mod- 
els which match the observed evolution of the 
more powerful sources. 

3.1. Simulating the radio source sky 

Simulations of the extragalactic radio sky, 
based on number density predictions from mod- 
els of the evolution of the radio luminosity func- 
tion allow us to predict typical source distribu- 
tions and natural confusion limit to projected 
SKA flux density limits. We show some exam- 
ples here for a lambda-dominated 'WMAP' cos- 
mology (f2 m =0.23, ^a= 0.73 and Hq= 71 km s _1 
Mpc" 1 ). 

We start by assuming that the radio sky com- 
prises three populations of sources, namely FRI, 
FRII and star-forming galaxies, and the LRLF 
and evolution for each can be determined: For 
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the star-forming galaxy population we adopt the 
local radio luminosity function from the 2dFGRS- 
NVSS galaxy sample at 1.4 GHz We adopt 
parameterised number-density luminosity evolu- 
tion for the star-forming population as deter- 
mined from the HDF-N and SSA13 fields Q3|. 

The local radio luminosity functions and evo- 
lution for the FRI and FRII radio galaxies are 
derived using the methodology of Jackson & 
Wall (1999). In summary, we fit exponential 
luminosity-dependent density evolution (LDDE) 
to the observed source count at 151 MHz. We 
find a best-fit to the observed 151 MHz source 
count for LDDE of the FRII population, cf) z = 
F(P,z).<fio, where 4>q is the LRLF. For the 
function F(P,z) we fit A/ maa; =15.70, z c =7.425, 
Pi=25.94 P 2 =27.94 for the FRII population, cou- 
pled with mild (M = —3.06) density evolution of 
the FRI population. 

From the three evolving RLFs we can produce 
simulated sky regions at a range of SKA frequen- 
cies: At 1.4 GHz we use the evolutionary fit at 
151 MHz plus 'beaming factors' to account for 
the quasar and BL Lac populations (Fig.l). For a 
single sky simulation we randomly position each 
source within a specified sky area. Each source 
is then randomly oriented, both on and into the 
plane. FRI and FRII radio galaxies are randomly 
assigned a total intrinsic size (lobe-to-lobe) be- 
tween 50 kpc and 800 kpc (as observed for the 
powerful 3CR sources, JE|)> assuming no source 
size evolution with redshift. FRIs and FRIIs 
are modeled as double-lobe structures, with each 
lobe being half the intrinsic total size adjusted for 
source orientation and redshift, i.e. foreshortened 
and flattened. Doppler-beamed versions of the 
FRI and FRIIs are set to be point sources. Star- 
forming galaxies are modelled as a circular disks 
with a randomly assigned intrinsic disk diameter 
between 10 and 100 kpc: these sizes are generous 
given the increasing evidence that spiral disks 
at high redshift are significantly smaller (e.g. 
|12|1. These disks are oriented randomly on the 
plane of the sky. Examples are shown here (Fig. 
2), a wider range of simulated sky images, at 
151 MHz, 325 MHz an d 1.4 GHz, are available at 
http: / /www. atnf.csiro.au/ ~cjackson/lowfreq-sky 

We can predict the number of sources per 



square degree and the fraction of sky covered by 
each population at 1.4 GHz as shown in Tabled 
We estimate the minimum resolution which would 
give an average inter-source spacing of 10 pri- 
mary beams on the sky, allowing for the vary- 
ing set of source sizes in a single sky simulation. 
Between 51.4 GHz= 10 /xJy and 100 nJy the mini- 
mum resolution requirement to avoid natural con- 
fusion ranges from 0.86 arc-seconds to 0.07 arc- 
seconds. Also of interest for deep SKA observa- 
tions is the natural confusion rate, i.e. the pro- 
portion of sources which overlap each other on 
the plane of the sky. We find that this varies from 
~10% of all sources at 10 /iJy, rising to ~60% at 
100 nJy. This means that at a flux density limit 
of 5i.4 ghz= 100 nJy, more than half of all de- 
tected sources will be in line-of-sight coincidence 
with another source. 

Table shows the source distributions by red- 
shifts: here we see that both the FRI and FRII 
populations extend to very high redshift (z ^8) at 
Si. 4 GHz < 10 /iJy. These sources would provide 
the foreground sources to HI line-of-sight exper- 
iments (Kanekar & Briggs, this volume and also 
®) 

4. SKA SPECIFICATION FOR AGN 
CONTINUUM SURVEY SCIENCE 

Based on the SKA sky predictions presented 
here, if the source count continues unchanged 
down to 100 nJy, there will be about 40 radio 
sources per square arcmin. We have modelled 
these sources with finite extents, such that these 
faint radio sources have dimensions comparable 
with those of the optical discs of galaxies. At the 
full sensitivity of the SKA, individual sources will 
appear blended at the faintest flux density levels. 

Resolutions of at least 0.05 arcsec at 1.4 GHz 
(i.e. 1000 km baselines) are necessary to sepa- 
rate individual sources with a baseline distribu- 
tion as specified in the SKA Science Requirements 
specification 15 . However, to image sources 
with sufficient detail to determine if the emis- 
sion comes from the disk, star formation or an 
AGN will require higher resolutions. It is also 
critical that the SKA resolution at least matches 
that of complementary facilities (e.g. JWST) for 
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Figure 1. Normalised observed and model differential source counts at f.4 GHz. The model count 
comprises 3 populations: FRlf radio galaxies (dashed line), FR1 radio galaxies (dotted) and star- forming 
galaxies (dot-dash). Data points are from: * Bridle et al. (1972), + Fomalont et al. (1974), A White 
et al. (1997), * Mitchell & Condon (1985), o Prandoni et al. (2000), A Gruppioni et al. (1997), and * 
Richards (2000). The polygon is the count estimate from the P(D) (background-deflection) analysis of 
Wall & Cooke (1975). 
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Table 1 

Predicted source surface densities at 1.4 GHz. 



Flux density limit 
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Table 2 

Predicted source distribution per square degree at 1.4 GHz. 



Flux density limit 
10 ^Jy 1 /xJy 100 nJy 
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Figure 2. Simulated sky images at 1.4 GHz, showing a 10 arcmin square region for flux density limits of 
1 /zJy (upper) and 100 nJy (lower). The difference between the pair of images at the same flux density 
illustrates cosmic variance, with nearby large-lobed sources dominating the right-hand plots. There are 
5 populations of sources shown: FRI galaxies (light, double-lobed), FRII galaxies (dark, double-lobed) , 
starforming galaxies (dark, single disk), beamed FRIIs (dark, star), and FRIs (light, star). 
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multi-wavelength studies. Also, detailed follow- 
up of SKA continuum radio surveys will be nec- 
essary for 'high precision' cosmology: Planck and 
its successor missions will require detailed fore- 
ground excision, where radio-detected AGN are 
contaminants to the CMB signal. 

Finally, it is vital that the SKA is capable 
of making simultaneous wide-band observations. 
These will be used to monitor the time- variability 
of radio sources, furthering the study of the mech- 
anisms driving intrinsic as well as external vari- 
ability, e.g. the intra-day variability (IDV) phe- 
nomenon. 
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